A&A manuscript no. 

(will be inserted by hand later ] 



Your thesaurus codes are: 

06 (03.11.1; 16.06.1; 19.06.1; 19.37.1; 19.53.1; 19.63.1) 



ASTRONOMY 

AND 
ASTROPHYSICS 



A Catalogue of Galaxies in the HDF-South: Photometry 
and Structural Parameters* 



o 
o 
o 

(N 
c3 



> 

o 

(N 

in 
o 
o 
o 



Of 
I 1 

O 



X 



M. Volonteri 2 , P. Saracco 1 , G. Chincarini 1,2 

1 Osservatorio Astronomico di Brera-Merate, Italy, 2 Dipartimento di Fisica, Univ. di Milano, Italy 
Received , 2000; accepted 



Abstract. We describe the construction of a catalogue 
of galaxies in the optical field of the Hubble Deep Field 
South. The HDF-S observations produced WFPC2 images 
in U, B, V, and I, the version 1 data have been made pub- 
lic on 23 November 1999. The effective field of view is 4.38 
arcmin 2 , and the 5a limiting magnitudes (in a FWHM 
aperture) are 28.87, 29.71, 30.19, 29.58 in the U, B, V 
and I bands, respectively. We created a catalogue for each 
pass-band (Isi4, V606> B450, U300), using simulations to ac- 
count for incompleteness and spurious sources contamina- 
tion. Along with photometry in all bands, we determined 
on the 1.814-selected catalogue (Isi4 < 26) structural pa- 
rameters, such as a metric apparent size, derived from the 
petrosian radius, an asymmetry index, light concentration 
indexes and the mean surface brightness within the pet- 
rosian radius. 

Key words: cosmology: observations - galaxies: photom- 
etry - galaxies: fundamental parameters 



1. Introduction 

The Hubble Deep Field South (HDF-S) was observed in 
October 1998 by the Hubble Space Telescope (HST). It is 
the southern counterpart of the Hubble Deep Field North 
(HDF-N) and shares its characteristics of depth and spa- 
tial resolution. The HDF-S is a four arcmin 2 survey lo- 
cated at RA 22h 32m 56s, DEC -60° 33' 02", observed 
during 150 orbits with the Wide Field Planetary Camera 

2 (WFPC2). The WFPC2 detector is composed by 4 chips: 
the Planetary Camera (PC) with higher resolution (0.05 
arcsec/pixel) and smaller field of view (35x35 arcsec 2 ) and 

3 Wide Field Cameras (WF), with a spatial resolution of 
0.10 arcsec/pixel and a field of view of 77x77 arcsec 2 . We 
will refer to the area of sky observed by the Planetary 
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Camera as PC field, and to the one observed by the three 
Wide Field Cameras as WF field. 

HDF-S images cover a wavelength range from the ul- 
traviolet to near-infrared (4 broad-band filters roughly 
corresponding to standard UBVI). Images were taken in 
four filters: F300W, F450W, F606W, and F814W, by a 
dithering technique. Released images are a combination of 
all individual exposures, weighted for the background sig- 
nal, and resampled to a pixel scale of 0.0398 arcsec/pixel 
using the "Drizzle" package (Fruchter, Hook et al. 1997). 
Details about observations and data reduction may be 
found in the Hubble Deep Field South web page (http:// 
www.stsci.edu/ftp/science/hdfsouth/hdfs.html). 

The main goal of the two HDF campaigns is to study 
the characteristics of faint galaxies and provide constraints 
on models explaining the formation and evolution of galax- 
ies, in particular the excess of faint blue sources. It is there- 
fore crucial to have a reliable catalogue, both for sources 
detection and photometry. In this paper we present a 
galaxy catalogue derived from the public version 1 images 
in fits format of the WFPC2 data, released on 23 Novem- 
ber 1999, along with documentation about data reduction 
and absolute calibration in the different bands. 

This catalogue gives for each galaxy, besides photo- 
metric information in all the optical bands, metric size, 
mean surface brightness, asymmetry index and light con- 
centration indexes, all fundamental for different cosmolog- 
ical and evolutionary tests. 

In Section 2 we describe the applied detection pro- 
cedure and in Section 3 the technique used to compute 
magnitudes, while in Section 4 we deal with the selec- 
tion criteria adopted to extract the sample. In Section 5 
we summarize the estimation of the sample completeness, 
while we discuss the possible oversampling of sub-galactic 
structures, such as HII regions, in high resolution imaging 
in Section 6. Finally in Section 7 and 8 we describe the 
technique used to recover galaxy colours and structural 
parameters respectively. 
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2. Sources Detection 

Object detection was performed across the area observed 
in optical bands (WFPC2 field) using the software for au- 
tomatic sources detection SExtractor (Bertin and Arnouts 
1996). The software firsts smoothes the images and then 
applies a threshold to identify peaks. As a smoothing func- 
tion we have chosen a Gaussian function having FWHM 
equal to the one measured on the images (i.e. 0.16 arc- 
sec) and a detection threshold of 1 sigma per pixel, with 
a minimum detection area equal to the seeing disk, have 
been adopted to pick up objects. The software allows to 
examine every identified group of contiguous pixels to de- 
blend overlapping sources: two sources belonging to the 
same group are considered different objects if they differ 
of 5 magnitudes or less (deblend-mincont=0. 01, deblend- 
ntresh=32). 

This choice of low thresholds aimed at obtaining a raw 
catalogue not biased against very faint or small source, in 
fact the first sample consisted of 6093, 4747, 9850, 5229 de- 
tections in the F300W, F450W, F606W and F814W band 
respectively. From this sample we then extracted a cata- 
logue optimizing completeness while minimizing the num- 
ber of spurious detections by selection criteria based on 
the signal-to-noise ratio of sources and simulations (as ex- 
plained in Section 4). 

3. Photometry 

Galaxies in this survey span a wide range both in redshift 
and in physical size. A fixed aperture would measure dif- 
ferent fractions of flux for galaxies dissimilar in shape and 
size and at different redshifts, while isophotal magnitudes 
suffer from the (1 + z) 5 cosmological dimming in surface 
brightness. 

"Pseudo-total" magnitudes were therefore estimated 
using the method of Djorgovski ct al. (1995) and Smail 
et al. (1995): they assigned the isophotal magnitude to 
sources with isophotal diameter larger than 9\ ss 2 — 3 
FWHM, while smaller sources are assigned an aperture 
corrected magnitude, that is the magnitude within an 
aperture Q\ corrected to the magnitude within 6* 2 , larger 
than 6\\ m = m(91) + Am, Am =< m(92) - m(91) >. 
In literature the choice of 9\ and #2 is based on multi- 
ples of the FWHM of images. In the HDF-S the excellent 
seeing needs a different approach, we therefore estimated 
magnitudes for our sample on the basis of the following 
steps: 

— For large sources, i.e. those sources having an isopho- 
tal diameter D iso > 61, and for "blended" sources, 
as flagged by SExtractor, we selected the SExtractor 
isophotal corrected magnitude; 

— to small sources, having D iso < 9\, we assigned an 
aperture corrected magnitude (estimated within 9\ 
and then corrected by Am to 2 , being 9\ < 9 2 ) 



Table 1. Residuals nikr-mdj vs. uik r for the Vgo6 selected 
sample, if 9\ = 1.20 arcsec and #2 defined as a function of 



Filter 


62 (arcsec) 


Am 


med(m fer -m dj ) 


X 


F300W 


1.59 


0.132 


0.08±0.40 


1.01 


F450W 


1.85 


0.163 


0.11±0.15 


1.04 


F606W 


2.15 


0.080 


0.06±0.10 


0.99 


F814W 


2.01 


0.145 


0.12±0.12 


1.05 



B\ has been defined as the minimum apparent diame- 
ter of a galaxy having an effective diameter r e = 10 Kpc. 
Hereafter we use a A = cosmology, with qo — 0.5 and 
Ho = 50 kms _1 Mpc _1 unless differently specified. With 
this choice 9\ = 1.2 arcsec. Since the correction Am is 
measured on a subsamplc of relatively bright galaxies, 
we defined an area for each band, Ago, such that 90% 
of galaxies belonging to the djorg subsample has isopho- 
tal area smaller than ^90 • #2 is defined as the diameter 
corresponding to a circle of area Ag . 

At bright magnitudes Kron's technique, based on an 
adaptive aperture, r\ = '^r'i(r) ' §i ves very good results. 
Kron (1980) and Bertin & Arnouts (1996) demonstrated 
that a photometry within an adaptive aperture (2.5ri) 
is expected to measure a fraction of the total flux be- 
tween 0.9 and 0.94. We then chose to use Kron's magni- 
tude (mfe r ) as a reference in order to test our method: if 
the flux fraction measured by Kron's technique is 0.94, the 
fraction estimated by our "pseudo-total" magnitude (m<#) 
is x — 0.94 • lQ 0A ( m kr-m d j) _ Q ur procedure is intended to 
correct the systematic underestimate (6-10%) of total flux 
typical of Kron's technique, which may be important for 
very faint sources. 

Table 1 and Figure 1 show clearly that statistically 
Am corrects for the flux underestimate typical of Kron's 
magnitudes. Moreover our estimate of total magnitudes 
has a narrower distribution at low S/N than Kron mag- 
nitude, see Figure 3, and in a plot magnitude-isophotal 
area (Figure 2) it is not evident any discontinuity in the 
passage between large and small sources (i.e. sources with 
9 > #1 or vice versa). This test confirms the validity of 
our choice of 9\. 

If not differently specified, magnitudes are expressed 
in the AB system, that is a system based on a spectrum 
which is flat in f„: m=-2.51og/ I/ - 48.60 (Okc 1974) . 

4. Selection of the Catalogue 

Our low detection threshold led to a large number of detec- 
tions (more than 4500 on less than 5 arcmin 2 ): it is there- 
fore necessary to evaluate the number of spurious sources. 
The depth and coverage for each filter is not homogeneous 
in the field of view, due to the variety of pointings that 
were combined together. The image depth fades towards 
the edges of the area covered, as well as in a cross-shaped 
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Fig. 1. Residuals mfc r -m<g vs. nifc r for the V606 selected 
sample, if Q\ = 1.20 arcscc, 8 2 as a function of Ago. 



Fig. 3. Kron's magnitude and our pseudototal magnitude 
vs S/N. At low S/N pseudototal magnitude has a narrower 
distribution than Kron's magnitude. 




DC 26 

•flat' 

25 • ".••fTVA'i • 



200 400 600 800 1000 

1(7 isophotal area 

Fig. 2. Plot of pseudo-total magnitude vs isophotal area: 
it is not evident any discontinuity in the passage between 
large and small sources (i.e. sources with 6 > 6\, corre- 
sponding to isophotal area > 706 pixel 2 ). This confirms 
our choice of 9\ . 



area between detectors which received much lower cover- 
age than the central region of each chip. As a consequence 
of the decreased image quality, the outer regions of each 
image are less reliable: in fact sky RMS is higher than 
the average one and some sources near the edges revealed 
at inspection to be spurious. We selected reliable sources 



by means of criteria based on S/N ratio and comparisons 
with simulations. 

The " Drizzle" algorithm (Variable-Pixel Linear Recon- 
struction) used to combine the various pointings preserves 
photometry and resolution and removes the effects of ge- 
ometric distortion, but it causes adjacent pixels to be cor- 
related. The pixel-to-pixel noise (cr s fe y )therefore underes- 
timates the true noise of a larger area by a factor 1.9. 
The noise measured on PC field is greater by a factor 
w 2 than that measured on the WF area. The S/N is then 
computed by a semi-empirical model (Pozzetti ct al. 1998, 
Williams ct al. 1996): S/N=R/cr tot , where R are net counts 
and of ot = R/(Tt exp ) + 2 • 1.9 2 a 2 ky A obj for WF sources, 
a 2 ot = R/(Tt exp ) + 2 • 4 • l.9o- 2 ky A ohj for PC sources. In the 
above formulas a s k y is the pixel-to-pixel sky RMS, t exp is 
the exposure time, T is the gain expressed in electrons per 
ADU, A bj and A s k v are respectively the object and the 
sky isophotal areas (in pixel) used to estimate the local 
background. 

The former term in the sum represents the Poissonian 
noise due to the source, the latter estimates statistical fluc- 
tuations in the mean value of sky, in the Poissonian ap- 
proximation. The factor 2 is linked to uncertainties in the 
determination of local background: the correct term would 
be 1. 9 2 T 2 a 2 ky A 2 b j/ A s k y , but since A s k v differs less than 
30% by the mean value of A ob j , we considered A s k y ss A Q bj ■ 

Our detection threshold corresponds to a minimum 
signal-to-noise ratio of S/Nwf = 1-34 and S/Npc = 0.67 
for the faintest sources detectable on the WF area and on 
the PC area respectively. 

In Table 2 we report for each filter the zeropoint (AB 
magnitude, Oke 1974), the sky RMS estimated by SEx- 
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Table 2. RMS sky values and 5a magnitude limit for the 
different bands. 



Filter 


zeropoint 


RMS 


nlHm 






(ADU/pix)xlO" 5 




F300W 


20.77 


1.674 


28.87 


F450W 


21.94 


2.284 


29.71 


F606W 


23.04 


4.126 


30.16 


F814W 


22.09 


2.960 


29.58 



tractor and the corresponding ha magnitude limit for a 
point source. 

We treated this problem statistically, in the hypothesis 
that noise is symmetrical with respect to the mean sky 
value. Operationally we have first created for each filter 
a noise frame by reversing the original images, in order 
to reveal the negative fluctuations and to make negative 
(i.e. undetectable) real sources (Saracco et al. 1999). Then 
we run SExtractor with the same detection parameter set 
used to search for sources in the original images detecting, 
by definition, only spurious sources. Applying a S/N=5 
cut off, after removing the edges of the images, we were 
able to reduce the spurious contamination to a negligible 
fraction (4%) on the WF area, while such a cut off is not 
able to reduce spurious detections to a reasonable level 
on the PC area being them more than 35%. In Figure 4 
the magnitude distribution of spurious sources obtained 
on the WF area and the PC area in the F606W band are 
shown. It is clear that the influence of spurious sources 
on the PC field is still remarkable after applying selection 
criteria, while the contamination is suppressed in the WF 
field. 

Thus, to avoid introducing such a large number of 
spurious by the PC data, we restricted the selection of 
sources to the central WF area only corresponding to 4.38 
arcmin 2 . On this area 450, 1153, 1694 and 1416 sources 
have been selected accordingly to the above criteria in the 
F300W, F450W, F606W and F814W band respectively, 
while the raw catalogues had 6093,4747, 9850, 5229 detec- 
tions in the same bands. 

In every magnitude bin we compared the number of 
sources in our final catalogue with the number of spuri- 
ous detections in order to get the contamination of false 
detection, shown in Table 3. 

We then removed stars from the sample by using the 
SExtractor morphological classifier. We defined as stars 
those sources brighter than Igi4=22 and having a value of 
the "stellarity" index larger than 0.9. This choice tends 
to underestimate stars both at faint magnitudes where 
no classification is considered, and at bright magnitudes 
where some fuzzy stars could be misclassified as galaxy. 
On the other hand this will ensure that our galaxy sample 
is not biased against compact galaxies. The star "clean- 
ing" procedure has classified and removed 14 stars at 
I8i4 < 22 in agreement with the number of stars found 



Table 3. Percentage of spurious detections in every mag- 
nitude bin, as estimated comparing detections on the on 
reversed (i.e. multiplied by -1) images with detections on 
real frames, that is percentage ratio of objects detected on 
reversed frames to sources detected in original frames in 
every magnitude bin. 



Pass-band 


26.25 


26.75 


27.25 


27.75 


28.25 


28.75 


U300 


1 


5.8 










B450 








0.9 


1.9 


2.9 




V606 











0.3 


1 


3 


1814 





0.5 


0.8 


4.2 







800 




mag 

Fig. 4. Spurious sources, defined as detections on the re- 
versed (i.e. multiplied by -1) V606 image. The upper panel 
refers to the WF field, the lower panel to the PC field. 
The solid line represents the initial detections, the dashed 
line represents spurious sources left after applying our se- 
lection criteria. 

in the HDF-N by Mendez et al. (1998) to this depth and 
in excess by a factor of two with respect to the prediction 
of the galaxy model of Bahcall & Soneira (1981). 

5. Completeness Correction 

Completeness correction for faint undetected sources 
strongly depends on the source apparent spatial structure 
besides to their magnitude. To estimate completeness via 
simulations we must account for sub-galactic structures. 
These could be quite common at high rcdshift and de- 
tectable by HST, and for the morphology which is that of 
the co-moving UV and B pass-bands and hence is strongly 
affected by star formation episodes. These features im- 
ply that "typical" profiles of galaxies are not able to well 
describe the shapes of a lot of galaxies in the HDF-S. 
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Thus, in order to reproduce the manifold of shapes which 
characterizes sources in the HDF-S we generated a set of 
simulated frames by directly dimming the original frames 
themselves by various factors while keeping constant the 
RMS (Saracco et al., 2000). This procedure has allowed 
us to avoid any assumption on the source profile while 
providing an artificial fair dimmed sample in a real back- 
ground noise. We thus define the correction factor c as the 
mean number of dimmed galaxies which should enter the 
fainter magnitude bin over the mean number of detected 
ones. It represents the inverse of the fraction of galaxies 
undetected in each bin. If m is the number of galaxies de- 
tected in the i-th bin of the original catalogue and TOj_|_i 
is the number of sources in the i+l-th bin of the simu- 
lated catalogue, the correction factor Cj+i corresponding 
to the ratio between the expected number of galaxies (rii) 
and the number of galaxies recovered ( rrii+i). The "true" 
number of galaxies in the i+l-th bin of the original cata- 
logue is then Ni + i= rii+i • 5j+i, where c is the mean over 
different simulations for the same frame. 

When dimming fluxes by a factor F=10 OAy , magni- 
tudes are y mag fainter, but also noise is lowered: if a is 
the sky RMS on original images, the dimmed frames have 
a RMS aia/F. We then added a frame of pure Poissonian 
noise, with a sky RMS 



2 _ 2 

noise ® 



1 



(a/F)W(l-- 



The final images have, by construction, the correct RMS. 
By choosing y=0.5 mag, the artificial noise added is 
o'noise = 0.6cr, i.e. about one third of the final image is 
due to pure Poissonian noise. As showed in Figures 5-6 
and in Table 4, this choice allowed us to correct up to 
V606=29, with brighter limits in the other bands. 

We tried also y—1 mag, but the estimated incomplete- 
ness was catastrophic, corresponding to 96%. This result 
seemed caused by the high simulated noise ( <7„ O i se =0.84 
cr, i.e. almost one half of noise is artificial), which caused 
also brighter bins to be incomplete. 

The completeness correction allowed the computation 
of differential number counts up to fainter magnitudes: we 
thus determined the slopes of the number counts relation. 
Our best fit gives j v = 0.47 ± 0.05, j B ~ 0.35 ± 0.02, 
7v - 0.28 ±0.01 and 7/ - 0.28 ±0.01 (see Volonteri et al. 
2000, for a detailed discussion). 

6. Galaxies or HII Regions? 

Colley et al. (1996) suggested that galaxies in the HDF-N 
may suffer from a wrong selection. High redshift galax- 
ies on optical images have a clumpy appearance: first the 
redshift moves the ultraviolet rest-frame light into the op- 
tical, so galaxies are observed in UV rest-frame where star- 
forming regions are more prominent, second the fraction 
of irregular galaxies is higher than locally (van den Bergh 
et al. 1996, Abraham et al. 1996) and a large number of 
galaxies display asymmetry and multiple structure. 



Table 4. Correction factor c and error, the factor c ac- 
counts for the galaxies missed in the detection due to the 
influence of noise. The "true" (corrected) number of galax- 
ies in each bin is N = cn, where n is the raw number of 
detections. 



Filter 


mag 


c 


0"e 


F300W 


26.75 


2.12 


0.07 


F450W 


27.25 


1.11 


0.05 




27.75 


1.56 


0.04 




28.25 


3.25 


0.11 


F606W 


28.25 


1.34 


0.02 




28.75 


2.67 


0.1 


F814W 


27.25 


1.09 


0.04 




27.75 


1.81 


0.11 


F110W 


27.25 


1.24 


0.14 




27.75 


2.56 


0.19 


F160W 


27.25 


1.93 


0.42 


F222M 


24.25 


4.09 


1.23 




22 23 24 25 26 27 28 
mag 

Fig. 5. V606 an d Isi4 number counts: raw counts are 
shown with empty symbols, corrected counts with filled 
symbols. 



Colley et al. (1996) stressed that compact high-redshift 
objects may appear more prominently than diffuse objects 
if their angular size is smaller than the point-spread func- 
tion (PSF). The cosmological dimming in surface bright- 
ness for these sources is less significant, leading to an en- 
hancement of compact sources over diffuse, resolved ob- 
jects. Since HST images have an excellent seeing, these 
clumps are not smoothed, so they may confuse detection 
algorithms (SExtractor, D AO FIND, FOCAS, etc), and be 
counted as several distinct faint sources. 

As suggested by Colley et al. (1996) in this case strong 
correlations between sources on scales < 10 kpc should 
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22 23 24 25 26 27 
mag 

Fig. 6. B450 and U300 number counts: raw counts are 
shown with empty symbols, corrected counts with filled 
symbols. 



be found. A good test is therefore the two-point angular 
correlation function. 

In order to check our sample against this effect, we 
computed this function on small angular scales: HII re- 
gions physical sizes are about 0.5 kpc, so a wrong selec- 
tion of the catalogue should bring a positive peak in the 
two-point angular correlation function at 0.25 — 1 arcsec. 
This scale corresponds to sizes less than 10 kpc for a wide 
range in redshift (0.8 < z < 3.5). 

We computed the correlation function by comparing 
the number of data pairs at given angular separation to 
the number of data-random simulated pairs at the same 
separation (Davis & Peebles, 1983). At fixed 9 the sum of 
data pairs is 



where Si e Sj are delta-functions on i-th and j-th galaxies 
positions. The sum over i is over all sources in the sample, 
and the sum over j includes only objects within a distance 
9 from particle i. 

We then created a random sample and computed the 
cross count sum 



where Si is as before, and <5j* is the delta-function for po- 
sitions of objects in the random sample within a distance 
9 from particle i. 



The resulting two-point angular correlation function is 
given by: 

n R DD{9) 

where — is the ratio of the mean density of random and 

no ■' 

data samples respectively. 

We analyzed our sample in I 8 i 4 and B 450 bands. The 
Isi4-band catalogue should suffer less from the effects de- 
scribed above, being selected in the reddest filter, the vice 
versa is true for the B45o-band catalogue. 

The two-point angular correlation function is compat- 
ible with zero in each band if computed on the whole cat- 
alogue. Following Colley et al. (1996) we selected "small 
galaxies" as defined by T> = ^/ab < 0.2 arcsec, where a, 
b are the intensity-weighted second moments. At 9 = 0.8 
arcsec it is evident a peak in the function, w(9) — 0.67 ± 
0.59 (Figure 7). This feature may not be significant due 
to large errors. 

However about 20-30% of sources in the B45o-band 
catalogue have separation< 1 arcsec. We therefore an- 
alyzed these sources, by cross-correlating Isi4 and B450 
catalogues. In the B4 50 -band catalogue we selected pairs 
with a separation < 1 arcsec which were not included in 
the Igi 4 -band catalogue. These objects were single sources 
splitted in the B45o-band (with B450 ~ 27-29), correspond- 
ing to a single detection in the Isi4-band. We then used 
SExtractor on the B450 frame, after choosing a higher 
deblend-mincont=0.1. 87 sources, with 21 < -B450 < 26, 
corresponding to about 7% of the whole sample, were then 
considered as single galaxies. 

This simple analysis shows that a remarkable fraction 
of sources in the HDF-S has a neighbour at small angu- 
lar distance. It should be pointed out that the HDF-S is 
a pencil-beam survey, and projection effects may be size- 
able. 



7. Optical Colours 

The estimate of galaxy colours could be biased towards 
the selection band and affected by an aperture effect. In 
order to test for this, we first determined colours simply 
by using SExtractor in the so-called double image mode on 
the l8i 4 -selected catalogue, we measured B 450 magnitudes 
within the la isophotal area determined in the Isi 4 band 
(B co ;). The (B-I)/ colour is then the difference between 
B co ; and the la isophotal magnitude in the Isi 4 band. 
We then measured the (B-T)b colour for the B-selected 
catalogue using the same method, but with the B450 band 
la isophotal area as a reference for both the Isi4 (Toz) and 
B450 magnitudes. 

For every source the B-I colours measured on the ba- 
sis of the B450 or Isi4 area are different. In Figure 8 we 
show the residuals (B-I)/-(B-I)s vs 1814. (B-I)s is slightly 
redder than (B-I)/), probably because the B45o-band 
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Fig. 7. The two-point angular correlation function w(9) 
estimated in the B450 band on small angular scales, for 
"small" sources, as denned by V = \fab < 0.2 arcsec, 
where a, b are the intensity-weighted second moments. 
The function is compatible with zero, but, within the error 
bars, it is possible to notice a peak at = 0.8 arcsec 
(corresponding to a physical distance of 4-6 kpc for 0.3< 
z <5). 

isophote underestimates the Isi4 band flux (broadly speak- 
ing Igi4 isophotal areas are larger than B450 isophotal ar- 
eas). 

In order to measure colours unbiased by the selection 
band we chose a different approach, that is we created a 
meta-image UBVI by summing all four frames, after nor- 
malizing each to have the same rms sky noise. We then 
run SExtractor in the double image mode: detection and 
isophote boundaries were measured on the combined im- 
age, while isophotal magnitudes were measured on U300, 
B450, V606j 1814 images individually. Using this proce- 
dure (Moustakas et al. 1997) both object detection and 
isophote determination are based on the summed image, 
and isophotcs are not biased towards any of the bands. 

We finally cross-correlated the catalogue obtained from 
the combined image with the Isi4 sample selected accord- 
ing to our criteria (see Image Analysis). We assigned a 
lower limit in magnitude to sources undetected in any of 
other bands (that is missing in our final U300, B450, V606 
catalogues). The limiting magnitude is the 5a isophotal 
magnitude within the isophote measured in the combined 
image. 

8. Structural Parameters 

For a subsample (Isi4< 26) we computed size, surface 
brightness, light concentration index, asymmetry index. 
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Fig. 8. Difference between (B-I) computed in Isi4 or B450 
bands vs Isi4. Colours are measured having as a refer- 
ence the lcr I 8 i4 and B450 isophotal areas respectively. 
(B-I)j-(B-I)b correlates with the Isi4 magnitude: (B-I)s 
is slightly redder than (B-I)/ at faint magnitudes, prob- 
ably because the B4 50 -band isophote underestimates the 
Isi4 band flux. 

The measure of the size of distant galaxies needs the 
use of a metric radius, that is a size independent on the 
redshift and light profile. We estimated a metric radius 
based on the Petrosian function (Petrosian, 1976), defined 
(Kron, 1995) as 

l dlnf(g) 

where 1(9) is the light growth curve. This function has the 
property: 

< i >e 

where 1(9) is the surface brightness at the radius 9 and 
< I >e is the mean surface brightness within 9. We de- 
fined for each galaxy the angular size as the value of #0.5 
such that 7/(0) =0.5 (Bershady et al., 1998, Saracco et al, 
1999). In order to determine the function 7/(0) we obtained 
the intensity profile, through multi-aperture photometry 
at equispaced (0.04 pixel) diameters, and subsequently in- 
terpolated it by a spline fit. There are sources whose n(9) 
is always larger than 0.5 (considering as a border the aper- 
ture i — th such that magi + i > magi); in these cases we 
defined 9 n , 5 = — 1, and also all quantities linked to O .5 
were arbitrarily set —1. We also measured the effective 
radius (half-light radius, r e //); for Isi4< 26 the relation 
9 vs r e ff is well fitted by O .5 = ^-^ r eff- After determin- 
ing O .5, we computed for each galaxy the mean surface 
brightness within 0q,5. 
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Abraham et al. (1994) and Abraham et al. (1996) 
showed that two indexes, namely an asymmetry index (A) 
and a central concentration index (C a h r ), are very useful in 
order to estimate a quantitative galaxy morphology. The 
former is determined by rotating the galaxy by 180° and 
subtracting the resulting image from the original one. The 
asymmetry index is given by the sum of absolute values of 
the pixels in the residual image, normalized by the sum of 
the absolute value of the pixels in the original image and 
corrected for the intrinsic asymmetry of the background. 
The concentration index is given by the ratio of fluxes in 
two isophotes, based on the analysis of light profiles. The 
measure of these indexes is independent of colour, though 
they correlate well with optical colours. 

We therefore computed both the asymmetry and the 
concentration indexes following Abraham et al. (1994) and 
Abraham et al. (1996), but also computed a different light 
concentration parameters. 

That is we computed 

r F(<8 . 5 ) 
° ? ' F(< 1.50o. 8 ) 

(Saracco et al., 1999), i.e. the ratio between the flux within 
radius #0.5 and the flux within 1.500.5- 

As discussed by Saracco et al. (1999), C v is indepen- 
dent from the redshift of the source, since it is related to 
a metric size and is independent of the asymptotic profile. 
Brinchmann et al. (1998) on the contrary pointed out that 
the central concentration C a br defined by Abraham et al. 
(1994, 1996) is redshift-dependent. C a br has been com- 
puted in order to classify our galaxies and compare our 
results with HDF-North. We tested the various parame- 
ters against apparent magnitude and checked any correla- 
tion with colours (Consclice et al. 2000). We selected two 
subsamples, composed respectively by galaxies with both 
B450-V6O6 and V606-I814 redder or bluer than a local ellip- 
tical or a local irregular galaxy. We used these subsamples 
as tracers of colours. 

The asymmetry index (Figure 9) seems not to be bi- 
ased: the faintest sources are on the whole more symmetric 
than brighter ones, but the presence of asymmetric objects 
also in the last bins suggest this feature to be linked to 
the nature of these galaxies. The trend towards high sym- 
metry may be due to the influence of noise, which makes 
the profile smoother. 

The central concentration index C a br defined by Abra- 
ham et al. (1994, 1996) seems to be biased against compact 
sources at faint magnitudes (Isu > 24.5), while C v does 
not correlate with apparent magnitude (Figures 10-11) 



9. Description of the Catalogue 

We presented a catalogue of galaxies in the HDF-S, cre- 
ated using the public version 1 images of the WFPC2 data. 
We created a catalogue for each pass-band (Isu, V606, 
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Fig. 9. Asymmetry index vs Isu magnitude: sources with 
both B450-V6O6 arid V606-Isi4 redder than local ellipticals 
are shown as stars, while sources bluer than local irregu- 
lars are shown as triangles. 
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Fig. 10. Central concentration index C v vs Isu magni- 
tude: sources with both B450-V6O6 and V606-Isi4 redder 
than local ellipticals are shown as stars, while sources 
bluer than local irregulars are shown as triangles. 

B450450, U300), which is available upon request. In V606> 
B450, U300 for each galaxy the catalogue reads: 

— isophotal flux and magnitude, with their errors (SEx- 
tractor) 

— isophotal corrected flux and magnitude, with their er- 
rors (SExtractor) 
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logue reads also the isophotal corrected flux and magni- 
tude, with their errors, a fixed circular aperture (0.6 arcsec 
radius) flux and magnitude, with their errors, Kron's tech- 
nique flux and magnitude, with their errors, and SExtrac- 
tor best choice for photometry (Kron's or isophotal cor- 



rected). The full catalogue is available at bttp:/ /cdswcb.u- 



strsbg.fr/Abstract.html and via 
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Fig. 11. Central concentration index C a br vs Isi4 mag- 
nitude: sources with both B450-V6O6 and V606-Isi4 red- 
der than local ellipticals are shown as stars, while sources 
bluer than local irregulars are shown as triangles. 

— fixed circular aperture (0.6 arcsec radius) flux and 
magnitude, with their errors (SExtractor) 

— Kron's technique flux and magnitude, with their errors 
(SExtractor) 

— SExtractor best choice for photometry (Kron's or 
isophotal corrected) 

— isophotal area (pix 2 ) (SExtractor) 

— coordinates in pixels of the barycentre, related to 
versionl frames (SExtractor) 

— right ascension and declination (SExtractor) 

— second order moments of the light distribution (SEx- 
tractor) 

— semi-major and semi-minor axis lenghts in pixels (SEx- 
tractor) 

— flags produced by detection and measurement pro- 
cesses (SExtractor) 

— stellar classification (SExtractor) 

— signal-to-noise ratio (Section 4) 

— pseudo-total magnitude and error (Section 3) 

For the Isi4-selected catalogue besides the above en- 
tries, we estimated also the colours (Section 7) and for 
Ism < 26, the petrosian radius (arcsec), the mean surface 
brightness within the petrosian radius (mag arcsec -2 ), 
light concentration indexes, that is C n and C a br, and the 
asymmetry index as computed by Abraham software (Sec- 
tion 8) . In Table 5 a small part of the catalogue is shown. 

Appendix A: The Catalogue 

Example of a part of the catalogue: we show only some 
photometric entries for sake of concision. The real cata- 
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